Critical Percolation and Transport in Nearly One Dimension 
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A random hopping on a fractal network with dimension slightly above one, d — 1 + e, is considered 
as a model of transport for conducting polymers with nonmetallic conductivity. Within the real space 
renormalization group method of Migdal and Kadanoff, the critical behavior near the percolation 
threshold is studied. In contrast to a conventional regular expansion in e, the critical indices of 
correlation length, u = + 0{e~^^'), and of conductivity, t ~ exp(— 1 — 1/e), are found to be 
nonanalytic functions of e as e — > 0. Distribution for conductivity of the critical cluster is obtained 
to be gaussian with the relative width ~ exp( — 1/e). In case of variable range hopping an "1-d 
Mott's law" exp [-(Tt/T)^/^] dependence was found for the DC conductivity. It is shown, that the 
same type of strong temperature dependence is valid for the dielectric constant and the frequency- 
dependent conductivity, in agreement with experimental data for poorly conducting polymers. 
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Conducting polymers represent a large class of new 
materials with a great variety of transport properties 
[01 . The room-temperature conductivity (ctat) of some of 
them attains the metallic value and the temperature and 
frequency dependencies of their conductivity are closely 
to be metallic. The nature of the metallic phase is 
presently a subject of intensive study One point is 
that the metallic state is provided by strong interchain 
coupling in these materials In polymers with mod- 
erate aRT (of the order of several hundreds S/cm) the 
conductivity, as a rule, decreases with decreasing tem- 
perature. Because this decay follows a power law in a 
large temperature interval, presumably these materials 
are near the metal-insulator transition which happens at 
the critical interchain coupling. Poorly conducting sam- 
ples with a JIT of the order or less 1 S/cm demonstrate a 
behavior that can be classified as a dielectric one Q: it 
is similar to that observed in amorphous semiconductors. 
Namely, DC conductivity is strongly dependent on tem- 
perature and its best fit is given by an "1-d Mott law": 
Sdc oc exp — (Tq/T)^/^. For a variable range hopping 
mechanism of transport the temperature dependence of 
conductivity is [|,|): E^c oc exp -(To/T)i/('*+i), where 
d > 1 is the system's dimension. This is not true, how- 
ever, in a Id case for which it should be the Arhenius 
dependence, S^c oc exp— Tq/T. 

Experimental measurements of the microwave conduc- 
tivity and the dielectric constant of these poorly con- 
ducting samples revealed that both are strongly de- 
pendent upon temperature too, most probably accord- 
ing to the same 1-d Mott's law. The theory of hopping 
transport predicts, however, only a very weak power tem- 
perature dependence for the frequency-dependent con- 
ductivity and the dielectric constant in two- and three- 
dimensional systems [||. 

In the present Letter we suggest an explanation for 
these peculiar features of conducting polymers that is 
based on the specific structure of the polymer's network. 



In the stretched polyacetylene this network is formed by 
coupled polymer chains oriented along some direction. 
Electronic micrographs shows that in these substances 
polymeric chains are organized into fibrils |^ , which may 
be distinctly seen to be subdivided into smaller ones 
|]lO| . In a non-fibrillar form of conducting polymers, like 
polyaniline. X-rays data reveal the existence of highly or- 
dered "crystallinity regions" with metallic properties [|| . 
Therefore the whole network in the stretched polyani- 
line may be thought of as constructed from long one- 
dimensional polymer chains randomly coupled by metal- 
lic crystallinity islands of various sizes. The volume frac- 
tion of latter ones can be small. 

We assume here that a polymer structure represents 
nearly one- dimensional fractal. That means a specific 
kind of polymer's chains organization, defined in the fol- 
lowing way: Choose a three-dimensional cube with the 
edge L. Chains, which are coupled within this cube, form 
a set of bundles disconnected from each other. If for large 
enough L the cross-section of the maximum bundle is 
proportional to L*^, where < e < 2, then we shall call 
the system d* = 1 -t- e -dimensional. Obviously e = 
for purely one-dimensional systems (sets of uncoupled 
chains). Note, that if one assume chains to be connected 
either with low concentration of uncorrelated interchain 
links, or with weak links (their resistivities being high 
compared to intrachain ones in our example), then we 
are dealing with a quasi-one-dimensional system 
which is three-dimensional according to our definition. 

The transport mechanism in conducting polymers in 
the localized phase is assumed to be the variable range 
hopping type (VRH) ||ll[. The regular method to treat 
VRH models is the effective medium approximation 
which gives wrong results in the nearly-ld case. For ex- 
ample, for the percolation model it gives a threshold con- 
centration of broken bonds Ct ~ e, while c* « exp (—1/e), 
as we shall see later. The results for the critical expo- 
nents are also wrong in this case. We choose the following 
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approach. We will first study the nearly-ld percolation 
system exactly. The VRH model is reduced to the perco- 
lation problem by constructing the effective percolation 
lattice |,§. 

Our first aim is the study of the critical behavior of 
the conductivity near the percolation threshold in a d- 
dimensional lattice, where d is close to the lower critical 
dimensionality, i.e. d = 1 + e, e <C 1. The real space 
renormalization group of Migdal and Kadanoff (RGMK) 
[^-0, that is exact at d = 1, expected to be the appro- 
priate tool if d is close to one. This method was applied 
to the percolation conductivity problem by Kirkpatrick 
[ psf , but the case of dimensionality close to 1 did not get 
any special attention. 

The RGMK method may be formulated as follows: the 
d-dimensional hypercubic lattice is replaced by a (m, n) 
hierarchical structure (see below) with rn = n'^^^, and 
n — > 1 afterwards. The (m, n) hierarchical pseudolattice 
may be constructed by the infinite repetition of two sub- 
sequent steps, as illustrated in Fig||: 

a) Formation of a n-length chain from n bonds; 

b) Formation of a m-bundle with m chains in a cross- 
section. 

Bundles obtained by this way are used as elementary 
bonds at the next stage, etc. To have a continuous RG 
transformation instead of discrete one, one should pro- 
ceed to infinitesimal transformation, setting n = I + 
z/ , m = 1 + £//; — > 0, where e = d — 1. At this step 
we introduce the continuous scale variable A — exp (I'l), 
where I is the order of bonds in the hierarchical structure. 
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dom resistivities or conductivities, respectively. These 
functions are related by the integral transformation: 



oo 

f ds 

P{a) = 1 - ^Ji{2V^)Q{s), 



(2) 



where Ji(s) is the Bessel function. For the hierarchical 
structures, we arrive after the transition to a continuous 
transformation, at the following equation: 



A 



dQ{s,X) 



da 

7^ 



Ji (2 Vi^)P(cr, A) In P{<7, A) = 0. (3) 



Finite concentration of broken bonds Cb may be con- 
sidered by introducing the boundary conditions: 



(5(s, A) \s=o= 1 - Cfc 



P(fT,A) U=+oo-Cfc. (4) 



Setting in Eq. (|) s — > 0, we obtain the RGMK equations 
for broken bonds concentration, 



A 



dcb 
dA 



ecblncb - (1 - Cb)ln(l - Cft). 



(5) 



There are three fixed points, two stable ones, Cb = 
(regular lattice) and C{, = 1 (completely broken lattice), 
and one unstable fixed point Cb ~ Ct, where the threshold 
concentration, Cj (0 < Ct < 1) is given by 



ectlnct = (1 - Ct)ln(l - Ct). 



(6) 



At arbitrary values of e one can solve Eq. (|3|) only 
numerically. Numerical results concerning this solution 
under the threshold boundary conditions: Q \s=o= ^ — Ct, 
Q |s=+oo= 0, support the following assertion: the func- 
tion Q{s, A) eventually evolves at sufficiently large A into 
the automodeling form: 



Q{s,X) = Q{sX''), P(a,A) -P(aA-'^), 



(7) 



where a is the e-dependent exponent. This form ensures, 
together with Eq. (|^) , the power-law critical behavior of 
the conductivity E: 



FIG. 1. The two initial stages of (2,3) hierarchical struc- 
ture formation. 

It appears to be convenient to work with conductivity 
S and resistivity R = l/S distribution functions in the 
Laplace representation: 



P(cr) EE (exp(-aE)); Q{s) = (exp(-si?)). 



(1) 



This is justified by the fact that for the chains or the bun- 
dles formation we have simply sums of independent ran- 



E(r) 



Ct 



t — av. 



Ct 



(8) 



Here u is the critical exponent of the correlation length 
(the characteristic size of a finite connected cluster) that 
is determined from Eq. (^ by linearizing near the thresh- 
old point: 



= e(ln Ci + 1) + 1 + ln(l - ct), 



(9) 



where Ct is given by Eq. (|^). The critical index for the 
conductivity, as a function of e can be found only nu- 
merically, and the result is shown in Fig. ||. 
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FIG. 2. The conductivity exponent t as a, function of 
e — d ~ 1. The insert shows the comparison of numerical 
(sohd curve) and analytical results. 

In the limit of small e the problem may be treated an- 
alytically. In this case the threshold concentration and 
the critical index of the correlation length are given by: 
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(10) 



One can see, that, though the main part of v^^{t) depen- 
dence is linear, the remaining dependence on e is essen- 
tially nonanalytic, and therefore no regular e-expansion 
in 1 -|- e dimensions is possible. 

Let us look for the automodeling solution at e <C 1. 
Using the fact that for e = this solution is simply 
exp(— s) (a = 0), we assume that 



Q{x) = (1 - Cf)exp[- 



(11) 



as well as that a ^ Ct, v ^ Ct- Plugging Eqs. (|ll|,|7D and 
( |lO| ) into Eq. (^), linearizing in v, and neglecting higher- 
order terms over ct, such as av, Cj , Ctw, etc., we arrive at 
the following equation: 



^v" ~^v' + v = e 



where x — e^/{l + e). 



1 r= — exp 



(12) 



The solution of Eq. (|12|) with t;(0) = t;'(0) = can be 
found by the substitution: u(^) = S,w{^). We have: 



w 



'(0 = ^ e"'/^[l 



j ^e-^/(i+')Jr(2V^)] 



-ea [l-(l+Oe-«]}. 



(13) 



To eliminate here the term, growing at ^ 3> 1 as exp ^ 
exp(a::/e), the following condition should be fulfilled: 



ea = Ct 



1 - 



-'7/(i+e)j^(20y) 



(14) 



This equation determines the critical exponent for the 
conductivity t = av (see Eq. (H)). By combining Eqs. 
(|l^) and (|l|) we get: 



■ exp 



{d-iy 



■ exp 



1 



(15) 



The coefficients Vn of a Tailor expansion of v(x) at 
X = represent the cumulants of the distribution func- 
tion for rcsitivities, e.g. V2 = ({R"^) — (R)^^ I {R}^ ■ 
They are obviously of the order of ct ■ That means that 
the resistivity distribution at length scale A is of an al- 
most gaussian 5-function, centered at (R) cx A" and with 
the width cx A° exp(— 1/e). 

Using the scaling arguments []l6| , ^ one can derive the 
low-frequency part of the AC conductivity E (w, r) near 
the percolation threshold, if the critical exponent of the 
DC conductivity is known. System's conductivity at low 
frequencies uj ^ Wo (Wo is the hopping rate of the retain- 
ing bonds) and near the percolation threshold, r ^ 1, is 
supposed to be represented in the scaling form 



S(c.,t) 



en. 



"aj\T\'Wo9 



kT 



_S_ 

Wo 



-u/t 



(16) 



where S — —iuj, a\\ is the longitudal lattice constant, and 
t 



s + 1 



s = 2v- p. 



(17) 



Here s and u represent the critical exponents for the di- 
electric constant near the percolation threshold and the 
frequency-dependent conductivity at the threshold, resp., 
and (3 is the exponent of the percolation order parameter 
(the probability of the site to belong to the infinite clus- 
ter). The scaling function g{x) has the following asymp- 
totic behavior: 

as l^l < 1 ; 

-(-..., as a; > 1 ; 
. , as a; < 0, |a;| 1. 




3 



The critical exponent (3 in nearly one-dimensional sys- 
tems may be obtained to be: 



P = = — exp — 



It may be set equally to zero in all the following expres- 
sions due to a sufficiently high order smallness of Ct, so 
that we have: s ~ 2v k, 2/ e. 

Now let us consider a nearly one-dimensional system 
with variable range hopping along the chains. We as- 
sume [H that the hopping rates depends on the inter- 
site distances r^j = jr^ — rj\ and on the energy difference 
£ij = — £j| via fij = Tijja -f £ij/2kT (a is the lo- 
calization radius) as: Wij = uJoe-xp{—2fij). Assuming 
site positions and energies near the Fermi level homoge- 
neously distributed with the density Np, we arrive at the 
"two-dimensional Poisson" distribution for fij, F{f) = 



exp 



if/foY 



, with /o = (To/T)l/^ kTo = (Npa) \ 

It is known that the behavior of the conductiv- 
ity with continuously distributed hopping rates may be 
explained, at least qualitatively, if the result for the per- 
colation lattice is known. For this purpose, we choose 
some trial value of the hopping rates Wc- Then, the ini- 
tial system is replaced by the percolation network where 
the hopping rates less than Wc are set to zero (broken 
bonds), and all others are set equal to Wc- The conduc- 
tivity of such a system is obviously less than the initial 
one. If the trial value of Wc is chosen so that the conduc- 
tivity of the percolation system is maximal, one can hope 
that this value gives a good estimate for the conductivity 
of the real system. 

Performing this procedure, e.g. at zero frequency, one 
can find, that the DC conductivity is determined by the 
hops with fij very close to its threshold value ft, which is 
determined by F{ft) — Cf. From Eq. ( [lO| ) it follows that 

ft = fo/V^ = {TijTf'^, where Ti = To/e = (eiVFa)"'. 



Taking in Eq. (|16|) rig = NpkT and an — l/ric, we have: 



e'Wt 



^DC 



Nf [kTy Nf [kry 



■ exp 



{Ti/T) 



1/2 



(18) 



At the lowest frequencies (hydrodynamical region, uj <C 
bJh) the conductivity can be close to its DC value and 
may be represented as a power series m. —iui: T,[uj) ^ 
"^DC (1 — i^/i^h — ■■■), with a very small value for ujh- 
As a result, the static dielectric constant 77 becomes very 
large and strongly temperature-dependent 



Here, the temperature dependencies of the conductivity 
and of the dielectric constant remain essentially of the "1- 
d Mott's" type. If the temperature is relatively high, i.e. 
T > T2 {T2 = e^Ti/4), this dependence would continu- 
ously turn into one that is typical for d — 1 (see, e.g. O]) 
at frequencies of the order of ujm, ln( uJm/Wt) ^ — 2/e. 
At lower temperatures, T <$:. T2, Eq. ( |2C| ) is valid 
only until uj < \n{Q^/Wt) - -(1/e) ln(T2/T). At 
ujm UJ UJM, where \n{ujM/uJm) e(ri/T)^/^, we 
have the following dependence: 



Nf {kry 



{-iuj) 



1/2 



e(Ti/T)^ 
In {—iuj/ujn 



2/e 



(21) 



77 ~ Wt/iUh ~ (ri/r)i/^exp(2/e2) 



(19) 



which changes into one-dimensional one (some power of 
temperature) at uj ^ lum- 

In sum, by taking into account the specific fractal 
structure of polymer's network we obtained the very spe- 
cific temperature dependence for the static conductivity 
(Eq.(p^), dielectric constan t (Eq.(|l^ )) and frequency de- 
pendent conductivity (Eqs.( po| , pl[ )). The source of this 
dependence is the following. The low dimensional ran- 
dom system can be separated in weakly coupled clusters 
within which carriers occur to be locked. With temper- 
ature the size of clusters exponentially increases, since 
more space becomes accessible for carriers due to ther- 
mal activation. As a result the dielectric constant and 
the conductivity exhibit strong temperature dependen- 
cies. There is also the very strong frequence dispersion 
of conductivity starting with extremely low frequencies. 
Such a type of strong temperature and frequency depen- 
dence is actually observed in conducting polymers wth 
localized carriers 

Let us stress again that a similar dependence can not 
be obtained within the standard 2d or 3d hopping mod- 
els. By assuming the strong Coulomb repulsion between 
carriers, the "Id Mott's law" could be reproduced 
but it does not lead to any strong dependence for the 
dielectric constant and for the frequency-dependent con- 
ductivity. The reason is the following. In contrast to 
the low-dimensional case the clusters in 2d and 3d sys- 
tems prove to be more effectively coupled. Therefore, the 
large polarization of clusters does not happen because of 
leakage of carriers. Thus, our results support strongly 
the idea that the conducting polymers represent a low- 
dimensional substance even in the dielectric phase. 
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At u; 3> w/i the frequency dependence of conductivity is 
given by: 



S(w) 



Nf (kry 



■ exp 



X 1/2 / . X e/2' 



T 



Wt 



(20) 
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